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Abstract. We present in this paper a homogeneous pho- 
tospheric abundance study, on the basis of the analysis of 
high resolution optical spectra, of six post-AGB objects 
displaying a 21/ini circumstellar dust feature in their IR 
spectrum. The F-G spectral type of the 21/im stars make 
that a large range of elements including a wide variety of 
s-process elements, can be studied by their atomic lines. 
The high C/0-ratios together with the large overabun- 
dance of s-process elements prove that the objects are de- 
scendants of field carbon stars. We discuss in detail the s- 
process abundance distribution displayed by these 21fiia 
stars and conclude that the 3rd dredge-up efficiency is 
closely related to the strength of the integrated neutron 
irradiation. The expected anti-correlation of the neutron 
irradiation with metallicity, on the other hand, contains a 
large intrinsic scatter. Finally we compare our results with 
other intrinsic and extrinsic s-process enriched objects and 
conclude that the post-AGB stars offer very useful com- 
plementary data to constrain the evolutionary models of 
AGB nucleosynthesis and dredge-up processes. 
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1. Introduction 

A sub-class of post-AGB objects display in their IR spec- 
trum an emission feature at 21/im, first identified on IRAS 
LRS spectra by Kwok et al. (1989). It was realised that the 
feature only occurs in spectra of C-rich post-AGB stars as 
evidenced by the presence of PAH-emission in the mid-IR 
spectra and often strong HCN measurements. The car- 
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rier(s) of the 21/im feature is (are) still unidentified. Its 
excitation and/or formation is limited to a short evolu- 
tionary phase since the feature is not observed in IR spec- 
tra of AGB stars nor in those of planetary nebulae (e.g. 
Kwok et al. 1999; Yolk et al. 1999). The carbon rich na- 
ture of the circumstellar environment of these objects was 
corroborated by the detection of C2, C3 and CN bands in 
the optical spectra (Hrivnak 1995). Detailed analysis of 
high-resolution spectra of these bands indicates that they 
are formed in a relatively narrow shell in the circumstellar 
envelope (Bakker et al. 1996; 1997). 

Several studies of the photospheric composition con- 
firmed the post-3rd dredge-up character of individual 
21/im stars: they do not only display a high photospheric 
C/0 ratio but also a large overabundance of s-process el- 
ements (Van Winckel 1999 and references therein). The 
colour-temperature of the IR-excess, the high luminosity 
and actual spectral type together with the chemical com- 
position of both the photosphere and circumstellar mate- 
rial are all observational evidence that the 21//m stars are 
post-carbon stars. 

Although in recent years the efforts and results on the 
theoretical modelling of the structural and chemical evo- 
lutionary AGB models are impressive (see e.g. Straniero 
et al. 1995; Frost & Lattanzio 1996; Herwig et al., 1997; 
Mowlavi et al. 1998; Mowlavi 1999; Langer et al. 1999) the 
calculations are not only dependent on fundamental pa- 
rameters like initial mass and metallicity, but also on much 
less constrained or even free parameters like mass-loss his- 
tory and geometry, details of the dredge-up phenomena 
and thermal pulse nucleosynthesis, numerical treatment of 
convection, interpulse phase at which the post-AGB phase 
starts, engulfed proton profile into the intershell, etc. 

There is now general agreement that the neutrons, 
needed for the s-process nucleosynthesis in AGB stars, arc 
mainly produced by the ^'^C(a,n)^^0 reaction with a pos- 
sible contribution from the minor activation of the reac- 
tion ^^Ne(Q:,n)^^Mg (e.g. Smith & Lambert 1986; Jorissen 
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a primary source of "'^'^C. This ^'^C pocket is naturally build 
up after engulfment of free protons into the ^^C-rich in- 
tershell. The engulfed proton profile is, however, still a 
free parameter in the predictions. Straniero et al. (1995) 
proposed that the s-process synthesis occurs during the 
intershell in a radiative environment and they computed 
that the neutron irradiation varies within the radiative 
intershell depending on the engulfed proton profile. 

The observational data guiding the theoretical chemi- 
cal evolutionary models come from the analysis of intrinsic 
AGB stars of the M-MS-S-SC-C star sequence, (e.g. Lam- 
bert et al. 1995 and references therein) which is thought 
to reflect, at least partly, the evolution on the AGB of sin- 
gle stars towards an increasing C/0 ratio. In addition, the 
enrichment in ^^C seems to correlate well with an increase 
in s-process abundances (e.g Smith & Lambert 1990). The 
study of the abundance distribution of different s-process 
elements is therefore a very useful tool to characterise the 
internal nucleosynthesis during and in between thermal 
pulses and to confine models of chemical evolution on the 
AGB. The strong molecular opacity and sometimes unsta- 
ble photospheres of AGB stars make a quantitative anal- 
ysis, however, difficult and some species even untraceable. 

A second group of objects with excess abundances of 
s-process elements is the Ba-stargroup. In the rest of the 
paper we use the Ba-stargroup as a generic term for all 
objects for which the excess is acquired by mass-transfer 
from the companion which is now a white dwarf. It in- 
cludes strong and mild Ba giants and Ba dwarfs, CH gi- 
ant and subgiants, extrinsic S stars and since recently also 
yellow symbiotics (e.g. review by Jorissen 1999). The char- 
acteristics of the s-process excesses in these stars refiect 
the internal chemical enrichment of the companion, and 
the abundances may have been diluted by mixing of un- 
processed material during the evolution of the gainer. The 
S-type stars in this group are normally labeled extrinsic 
S-stars and are characterised by the lack of Tc in their 
atmospheres (e.g. Jorissen et al. 1993) or their lower lu- 
minosity than intrinsic AGB S-stars (Van Eck et al. 1998). 

In this paper, we analyse another group of intrinsic s- 
process enriched objects, namely the 21/zm stars. With 
their F to G spectral types, a wide variety of chemi- 
cal species can be studied quantitatively by using photo- 
spheric atomic lines. This makes them ideal complemen- 
tary sources to constrain the 3rd dredge-up nucleosynthe- 
sis and dredge-up models. 

The chemical composition of individual 21/im objects 
is already reported in the literature (Klochkova 1995; Zacs 
et al. 1995; Van Winckel et al. 1996a; Reddy et al. 1997; 
Decin et al. 1998, Klochkova et al. 1999). Since the dif- 
ferent authors use different line-lists, spectral resolution 
and atomic data, it is often difficult to compare the re- 
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the different s-process elements and the integrated neu- 
tron irradiation. In Sect. 2 we give an overview on the 
data obtained and sketch the reduction method, while the 
details of the chemical analysis are given in Sect. 3. In 
Sect. 4 a synopsis of our abundance analysis is presented 
proving that the 21/im stars show clear chemical evidence 
for their post-AGB character. Sect. 5 is devoted to the 
discussion of the results of the individual objects and in 
Sect. 6 we focus on the heavy-element nucleosynthesis and 
the determination of the neutron exposure. Our findings 
on the post-AGB stars are then compared in Sect. 7 to 
the s-process abundances displayed in intrinsic AGB stars 
and extrinsic heavy-element enriched objects with a wide 
range in metallicity. We end with formulating the main 
conclusions in Sect. 8. 

2. Observations and reduction 

In Table |] we list the sample of 21/im stars discussed here 
together with the spectral domain observed and signal- 
to-noise obtained. Most spectra were obtained using the 
Utrecht Echelle Spectrograph (UES) mounted on the 4.2m 
William Hershel Telescope (WHT) on La Palma, Spain. 
We used the echelle with 31.6 lines/mm and the projected 
slit width was 1.1" on the sky yielding a resolution of 
around R = X/SX ^ 50 000) depending on the wavelength. 
Some of the UES spectra were kindly provided by Dr. 
Eric Bakker and used by him in his analysis of the opti- 
cal circumstellar molecular absorption bands (Bakker et 
al. 1996; 1997). For IRAS07134 the 3.5m New Technology 
Telescope (NTT) of the European Southern Observatory 
(ESO) was used in combination with the high-resolution 
echelle mode of EMMI using grating number 14. The pro- 
jected slit width was 1" on the sky, yielding a resolution 
of R ~ 60 000. Since most objects are heavily reddened, 
the S/N achieved depends strongly on wavelength and of 
course on the specific characteristics of the cross-dispersed 
spectrographs. The standard reduction included bias cor- 
rection, flat-fielding, cosmic hits identification and clean- 
ing and background subtraction. The extraction was done 
by a simple mean over the spatial profile and the wave- 
length calibration was performed by using Thorium- Argon 
lamp measurements. The different orders do not overlap 
from ~ 580 nm red-wards for the UES in combination with 
the 1024^ pixel TEK CCD and for 742 nm red-wards for 
EMMI with the 2048^ TEK CCD (ESO number 36). We 
normalised the spectra by dividing the individual orders 
by a smoothed spline function defined through interac- 
tively identified continuum points. We used the specific 
echelle context of MIDAS for the reduction. 

3. Analysis 
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Table 1. The programme stars: magnitude, coordinates, log of the observations and radial velocities. 



IRAS 


Other 
name 


Visual 
magnitude 
m(v) 


Equatorial 
coordinates 

O?2000 ^2000 


Galactic 
coordinates 
1 b 


04296+3429 




14.2 


04 32 56.6 


+34 36 11 


166.24 


-9.05 


05341+0852 




12.8 


05 36 54.2 


+08 54 10 


196.19 


-12.14 


07134+1005 


HD 56126 


8.3 


07 16 10.2 


+09 59 47 


206.75 


+9.99 


19500-1709 


HD 187885 


9.2 


19 52 52.7 


-17 01 50 


23.98 


-21.04 


22223+4327 




9.7 


22 24 30.7 


+43 43 03 


96.75 


-11.56 


23304+6147 




13.1 


23 32 45.0 


+62 03 49 


113.86 


+0.59 




Date and UT 


Telescope+ 
Spectrograph 


Sp. Range, 
(nm) 


S/N 


Whclio 

(kms^^) 


vlsk 
(kms^^) 


04296+3429 


22/2/1994 21:20 


WHT+UES 


556.5T022 


80 


-53 


-60 


05341+0852 


08/8/1995 05:15 


WHT+UES 


556.5T004 


50 


24 


8 


07134+1005 


29/9/1998 09:01 
24/2/1992 21:57 


NTT+EMMI 
WHT+UES 


398-662 
535T040 


blue : 160 

green : 210 

100 


86 

84 


71 
70 


19500-1709 


06/8/1995 23:55 
08/8/1995 00:17 
07/8/1995 01:34 


WHT+UES 
WHT+UES 
WHT+UES 


364-458 

453-680 

551-1004 


140 
240 
320 


12 
13 
12 


24 
25 
24 


22223+4327 


23/8/1994 01:56 
21/8/1994 22:45 


WHT+UES 
WHT+UES 


443.5-650.5 
556.5-1004 


180 
190 


-42 
-42 


-30 
-30 


23304+6147 


23/8/1994 03:48 
22/8/1994 01:09 


WHT+UES 
WHT+UES 


443.5-650.5 
556.5-1004 


50 
110 


-26 
-26 


-16 
-16 



line lists of the Vienna Atomic Line Database (VALD2, 
http://www. astro. univie.ac.at'^vald/, Kupka et al. 1999) 



were used, and estimated line strengths were computed. 
Next, the lines stronger than our detection limit were iden- 
tified in the measured spectra. 

The oscillator strengths (log((7/)-values) were taken 
from several sources: for Fe, we used the critical com- 
pilation of Lambert et al. (1996), completed with values 
of Blackweh et al. (1980) for Fell. Values for C, N and 
O were mainly taken from the Opacity Project (Hibbert 
et al. 1991, 1993; Biemont et al. 1991); values for the s- 
process elements from the Vienna Atomic Line Database 
(VALD2); values for the other elements mainly from The- 
venin (1989, 1990), but also from Wiese et al. (1966), 
Reader et al. (1980), Fuhr et al. (1988), Martin et al. 
(1988), Venn (1995) and Gonzalez et al. (1997). The hue 
list can be obtained from the authors upon request. 

3.2. Radial velocities 

The heliocentric radial velocities of the objects are given in 
Table n^ with an internal accuracy of ~1 km s~^ . Compar- 
ison with published velocities reveals no significant vari- 
ability in most objects (see Table |). For IRAS04296, an 
accurate velocity is hampered by our lower signal-to-noise. 
IRAS19500 and IRAS23304 are clearly not variable. For 
IRAS05341 and IRAS22223 we have not enough veloci- 
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itoring campaign, which spanned 440 days. They deduced 
radial velocity variations between 81.73 and 91.79 km s~^ 
for this pulsating star but did not find evidence for binary 
motion. The measurement of Parthasarathy et al. (1992) 
deviates significantly from this range but it is difficult to 
evaluate since they do not give the date of observation. 
From the reception date of the article itself, it could have 
been very well during the monitoring campaign of Lebre 
et al. 

We can conclude that for none of the stars there is 
evidence for binarity from the radial velocities obtained 
so far. 



3.3. Determination of atmospheric param,eters and 
abundances. 

In order to calculate abundances, we used the model at- 
mospheres of Kurucz (1993) on his CDROM nr. 13. These 
LTE atmospheres are uniquely determined by the effective 
temperature (Teff), gravity (log 5) and overall metallicity 
(Z). The models were computed using opacity distribution 
functions assuming a constant microturbulent velocity (^t) 
of 2 kms-i. 

The model parameters were based solely on the basis 
of the spectra, more specifically on the analysis of a set 
of well measured Fe I and Fe II lines with accurate atomic 
data. We used Fe for this purpose, not only because of 
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Table 2. Heliocentric radial velocities of the programme 
stars, chronologically tabulated. The data used for the de- 
termination of the velocities include CO, and CS rota- 
tional line emission and optical spectra (opt.). The veloc- 
ities we found are in bold. For these velocities, a detailed 
time indication can be found in Table n, together with 
their value with respect to the Local Standard of Rest 
(wlsr)- 



Table 3. Model parameters of the six programme stars. 



IRAS 


Whclio 

(kms"^) 


method 


date 

(mm/yy) 


ref. 


04296+3429 


-55 


CO 


12/88 


11 




-59 


CO 


01/89 


7 




-53 


opt. 


02/94 






-56 


opt. 


02/97 


3 


05341-^0852 


24 


opt. 


08/95 






25 


opt. 


12/96 


10 




26 


CO 


04/98 


12 


07134-H005 


87 


CO 


01/89 


7 




105 


opt. 


(?) 


9 




82-^92 


opt. 


01/91-04/92 


5 




86 


CO 


04/90-01/91 


1 




84" 


opt. 


02/92 






91 


opt. 


01/93 


2 




86 


opt. 


12/93&02/94 


8 




se*" 


opt. 


09/98 




19500-1709 


13 


CO 


03/86 


6 




13 


CO 


07/87 


4 




13 


CO 


07/90 


7 




12 


CO 


04/90-01/91 


1 




12 


opt. 


08/95 




22223+4327 


-37 


CO 


01/88 


6 




-42 


CO 


01/89 


7 




-42 


opt. 


08/94 




23304+6147 


-26 


CO 


12/88 


11 




-26 


CO 


01/88 


6 




-26 


CS 


01/89 


7 




-26 


opt. 


08/94 





(a) UES spectra (1992), (b) EMMI spectra (1998) 
ref. (1) Bujarrabal et al. 1992 (2) Klochkova 1995 (3) Klochkova 
et al. 1999 (4) Knapp et al. 1989 (5) Lebre et al. 1996 (6) Likkel 
et al. 1991 (7) Omont et al. 1993 (8) Oudmaijer & Bakker 
1994 (9) Parthasarathy et al. 1992 (10) Reddy et al. 1997 (11) 
Woodsworth et al. 1990 (12) Hrivnak & Kwok 1999 



face gravity log g by demanding ionization equilibrium for 
Fel and Fell; the microturbulent velocity £_t, by forcing 
the abundances to be independent of the reduced equiv- 
alent widths \og{Wx/X). Results of other authors or pho- 
tometric data were used as a first guess. The stars in the 
sample are ~1000K too cool to show Hcl-lines, so we 
could not obtain He- abundances, except for the hottest 
star (IRAS19500, Teff=8000K). We already reported the 
detection of a Hel-line (A587.565nm) for this star in Van 
Winckel et al. nflflfiaV The obta.ined nara.meters are listed 





Z 


Teff 

(K) 


logs 


(kms-i) 


IRAS04296 


-0.5 


7000 


1.0 


4.0 


IRAS05341 


-1.0 


6500 


1.0 


3.5 


IRAS07134 


-1.0 


7250 


0.5 


5.0 


IRAS 19500 


-0.5 


8000 


1.0 


6.0 


IRAS22223 


-0.5 


6500 


1.0 


5.5 


IRAS23304 


-1.0 


6750 


0.5 


3.0 
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Fig. 1. Sample spectra of the programme stars around 
the red carbon multiplet. The stars are velocity corrected. 
On top, the spectrum of the reference star HR 1865 is 
plotted. This massive supergiant (FOIb) has similar atmo- 
spheric parameters as the programme stars (Toff =7500 K, 
logg=2.0 and Ct=3.0 kms~\ see Decin et al. (1998)) but 
obviously no enrichment of helium burning products. 



and Van Winckel (1997) for an extensive error analysis and 
remind that the uncertainties, assuming LTE, are mainly 
induced by uncertain \og{gf) values and an uncertain tem- 
perature determination. Since the model atmospheres of 
Kurucz are listed with steps of ATcff = 250 K and Alog{g) 
= 0.5, that correspond to the uncertainties of the spectro- 
scopic parameter determinations, we did not interpolate 
between models. The possibility to infer an abundance of 
an ion critically depends on the presence of useful lines of 
that ion in the studied spectral range. We only used in this 
LTE analysis lines which appear unblended in our spectra 
and have an ennivalent w^idth smaller than 150 mA. 
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Fig. 2. Sample spectra of the programme stars and 
HR1865 around the oxygen triplet at A615.6nm. For the 
cooler stars, a reliable oxygen abundance is rather prob- 
lematic to derive, because of the small number of lines: 
only this multiplet is available for these stars. Hotter 
stars (IRAS07134 and IRAS19500) show more O hues. 
The equivalent width of the Ball-line exceeds our upper 
limit of 150 mA and can therefore not be used in a barium 
abundance calculation. A barium abundance is always dif- 
ficult to determine because this element shows only a few 
very strong saturated lines and weak optical lines with 
WA<150inA are very rare. 



very consistent C-abundances (^(CI) < 0.14, except for 
IRAS04296: cr(CI) = 0.20). Useful nitrogen lines (NI) are 
more difficult to obtain, but there are some weak lines 
in the red part of the spectrum; only for IRAS07134 we 
measured more than 5 weak lines. The oxygen abundances 
of the 21/im stars are in most cases more difficult to ob- 
tain. Since the O triplet at A777.4 nm is strongly non- 
LTE sensitive, weaker lines should be used. For the cooler 
stars in the sample, only the high excitation multiplet at 
615.6 nm can be used (Fig. 0), but since these lines are 
heavily blended, we used multiple Gaussian fitting to es- 
timate their equivalent widths. For the two hotter stars 
IRAS07134 and IRAS19500, there are more lines in the 
spectrum, so the 0-abundance is better established. 

From sodium (Na, Z=ll) to sulphur (S, Z=16), abun- 
dances are hard to obtain due to the lack of lines. Only 
two different sodium lines were found in the spectra 

{\EP,R9t^'^r^m anrl X.^^fiS S99 nm'l hnth nf tlipm VinAnntr 



dance calculation. Silicon (Si) has more lines, but it shows 
for some stars a rather large line-to-line scatter. 

For the heavier metals (Ca, Sc, Ti, Cr, Ni), more lines 
were found. Neutral calcium (Ca) is found in every star. 
Also scandium (Sc) and titanium (Ti) abundances were 
derived for almost every star, and this from lines of singly 
ionised atoms. Chromium (Cr) displays lines from neutral 
as well as singly ionised atoms, allowing an extra check 
for the ionization equilibrium (see Sect. 3.3). Nickel (Ni) 
is also found in every star, showing neutral lines in the 
cooler objects, ionised lines in the hotter ones. Other iron- 
peak elements found are vanadium (V), manganese (Mn) 
and zinc (Zn). 

For the 21^m stars which have F spectral types, the 
s-process elements are determined by lines from singly 
ionised atoms. The strontium (Sr) abundance is very dif- 
ficult to determine accurately since the resonance lines 
at A421.5nm and A407.7nm are heavily saturated and 
weaker optical lines are not present. Yttrium (Y) and zir- 
conium (Zr) abundances are more easily obtained. More- 
over, these elements show a large number of lines in the 
hotter stars of our sample. An accurate barium (Ba) abun- 
dance is, like for Sr, difficult to determine, since weak 
optical lines are absent and only a few very strong sat- 
urated lines are present (for example, see Fig. 0). Three 
Ball-Hnes were found with Wx < 150 mA: A416.600nm, 
A585.367nm and A873.776nm. Lanthanum (La), cerium 
(Ce) and neodymium (Nd) are easier to study and we 
obtained abundances for every star. Praseodymium (Pr) 
and samarium (Sm) have less lines and we were not able 
to determine abundances for these elements for every star. 
Also europium (Eu) lines are quite rare; we detected only 
three different EuII-lines in our spectra. Finally, we also 
detected a Hafnium (Hf, Z=72) hnc (A409.316nm) in the 
spectrum of IRAS07134 and IRAS19500 l^ig. ||). This line 
will be discussed in more detail in Sect. 

3.4- Hyperfine splitting 

We investigated the influence of hyperfine splitting (hfs) 
on some observed lines which are considered to be sensitive 
for this effect. Hfs causes extra desaturation but is only 
important for lines of elements with an odd atomic number 
and strong nucleon-electron interaction. 

Hfs will affect weak lines in their profile, but not in 
their equivalent width, leaving the derived abundance un- 
changed. For stronger lines on the flat part of the curve 
of growth, however, also the equivalent width will be af- 
fected. The hfs-influence critically depends on the temper- 
ature, since at higher temperatures the Doppler widths 
will wash out the effect of hfs on the line-profile. 

The Kurucz spectrum synthesis program "Synthe" 

^TCnrnr-z 1 QQ.'^'l ixrnc; nspH tn stnrlv ttip inflnpnrp nf hfs Wp 
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Fig. 3. The region around A679 nni for the six programme 
stars and HR1865. In this region, both hnes of hght 
s-process elements (Is) and of heavy s-process elements 
(hs) are present. Comparison with the normal supergiant 
HR 1865 proves that our stars are clearly s-process en- 
hanced. The ratio of the line strength of the Lall-line at 
A677.427nm and the Yll-line at A679.541 nm for each star 
gives a first very qualitative idea of the neutron exposure 
(which is discussed in Sect. ||). In particular, from this 
comparison we can expect a high neutron exposure for 
IRAS05341. 



ing lines which are considered to be among the most in- 
fluenced by hfs-splitting: one aluminium (A 394.401 nm), 
one scandium (A 660.460 nm) and two europium lines (A 
412.973 nm and A 664.506). The hsf-dccompositions were 
obtained from McWilham et al. (1995) and Biehl (1976). 

The influence of hfs on Wa of the synthetic lines 
was found to be noticeable, but spectrum synthesis made 
clear that the other sources of uncertainty (like continuum 
placement, errors on the atmospheric parameters etc.) 
are much more important in our programme stars. The 
small influence, even on the normally strongly affected Eu 
412.973 nm line may surprise, but this line is only observed 
in the hotter stars where it has an equivalent width of only 
32 mA at most. For the cooler 21/im objects this line was 
not observed but our calculations indicate that this strong 
line would be heavily affected indeed. 

Since line-profile fitting induces other uncertainties 
and since the effect on the weak lines used in our anal- 
ysis is very small, we did not include hfs profile fitting 



4. Synopsis of the abundance results 

In this section we discuss Table H and Fig. 0, which give 
an overview of the results of our abundance analysis. In 
the last column of Table El we list a typical abundance 
pattern for an unevolved metal deficient star (e.g. Lam- 
bert 1989; Wheeler et al. 1989; Edvardsson et al. 1993; 
McWilliam 1997) reflecting the chemical composition of 
our programme stars at the time of their formation. Com- 
paring these unprocessed abundances with the observed 
abundances permits us to detect possible changes in the 
chemical composition during the stars evolution. Applied 
to low-mass post-AGB stars, we expect the chemical sig- 
nature of the 3rd dredge-up without activation of hot- 
bottom burning i.e. the enrichment of the stellar photo- 
sphere with material exposed to He-burning during the 
TP-AGB phase. This signature includes: carbon and (to 
a lesser extent) oxygen enrichment, a C/0 number ratio 
higher than solar (C0/O0=O.51) and, above all, s-process 
enrichment. All 21^m sources studied here turn out to 
show almost all these enrichments, which make them a 
group of definite post 3rd dredge-up post-AGB stars. 



4-.1. Metallicity 

To calculate the metallicity, we used the iron abun- 
dance inferred from the Fell-lines since this is by far 
the dominant ionization level. The metallicity relative to 
the solar value ranges from -0.3 (IRAS22223) to -1.0 
(IRAS07134). This is a strong indication that our sam- 
ple indeed consists of an old and hence low-mass popu- 
lation. This interpretation is further strengthened by the 
high galactic latitudes of the stars, in the range |6|=9°- 
21°. An exception with respect to the latter criterion is 
IRAS23304, for which b = 0.6°. Its low-mass nature is, 
however, clear from the low metallicity ([Fe/H]=— 0.8). 



4-2. CNO abundance 

All stars display huge carbon-enrichments. This enrich- 
ment is beyond any doubt even when the internal scatter 
of the C-lines and the uncertainty induced by the model 
parameters are taken into account. As already pointed 
out in Sect. 3.2, the 0-abundance is in most cases less 



certain. Only for the two hotter objects IRAS07134 and 
IRAS19500 we feel confident that we detected a moder- 
ate but real 0-enrichment. The reliability of the oxygen 
abundances will be discussed in the following section, in 
which we treat the objects separately. 

Unfortunately, the errors on mainly the 0-abundance 
preclude accurate C/0 number ratio determination based 
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Table 4. Synopsis of the abundance results. The [a/Fe] value is the mean of the [el/Fe] values of the observed a- 
elemcnts Mg, Si, S, Ca and Ti; [Is/Fe] the mean of the [el/Fe] values of Sr, Y and Zr; [hs/Fe] the mean of [el/Fe] of 
Ba, La, Nd and Sm; [s/Fe] the mean of the former 7 s-process elements. 
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Jf..3. a-elements 

The [a/Fe] value is the mean of the [el/Fe] values of the 
observed a-elements. If we take into account a typical er- 
ror of 0.2 dex in the abundances, we may conclude that 
the a-overabundances reflect the chemical history of the 
Galaxy and are therefore not enhanced. 

The a-elements show, however, large differences in 
their [el/Fe] values within each star, up to an extreme 
value of 1 dex. Unfortunately, the lack of useful lines 
for the light a-elements prevent too strong conclusions 
whether this is intrinisic or due to non-LTE effects and/or 
unaccurate log{gf). The [a/Fe] value is therefore only an 
indicative number. 



4- 4- s-process enrichment 

The s-process elements observed in evolved stars can be di- 
vided into two groups: the light s-process elements around 
the magic neutron number 50 (Sr, Y, Zr) and the heavy 
s-process elements around the magic neutron number 82 
(Ba, La, Ce, Pr, Nd, Sm, (Eu)). In order to investigate the 
s-process abundances in more detail, four indexes are gen- 
erally defined: [s/Fe], [Is/Fe], [hs/Fe] and [hs/ls]. Which el- 
ements are taken into account to determine these indices 
is different from author to author and is mainly deter- 
mined by the possibility to compute accurate abundances 
of the different elements. To compare the results on 21/im 
stars with other s-process enriched objects, the same in- 
dex should be used. We follow the proposition of Busso et 
al. (1995) and define the Is-index as the mean of Sr, Y and 
Zr and the hs-index as the mean of Ba, La, Nd and Sm. 
Consequently, [s/Fe] is the mean of the 7 former elements 
and [hs/ls] = [hs/Fe] -[Is/Fe]. 

For unobserved species of the light s-process no correc- 
tion factor was taken into account since the odd-even effect 
is not very strong. For unobserved species of the hs index, 

fnr wrViirVi tVip nrlrl-p.vpn pflFprt is Trmr'li strnncrpr tVip nnrlp- 



This method was used to estimate the Sm-abundance for 
IRAS04296 and IRAS19500. Note that the [hs/ls] index 
is by its definition in principle independent of the total 
s-process enrichment ([s/Fe]). 

Regarding Table H, the overabundance of the s-process 
elements is beyond any doubt and is definitely the most 
convincing argument for the post 3rd dredge-up status of 
the studied stars. For unevolved objects in the same metal- 
licity range, one expects the s-process elements to scale 
with Fe: [s/Fe]=0 (see Wheeler et al. 1989). Here we find 
[s/Fe]>0.9 (IRAS22223), up to [s/Fe]=2.2 (IRAS05341). 

5. Individual objects 

The individual objects will be discussed in this section. In 
Table ^ we list the complete abundance analysis of each 
programme star. We want to stress that we will not repeat 
the conclusions already made in the previous section, but 
that we will focus on the peculiarities in the analysis of 
each star. 



5.1. IRAS 04296+3429 

The chemical analysis of this object has already been dis- 
cussed by our group (Decin et al. 1998). Nevertheless, 
there are some minor differences between the present anal- 
ysis and the previous one, as a consequence of an update 
of our line list. The main conclusions of Decin et al. (1998) 
remain valid. 



5.2. IRAS 05341+0852 

This object was first mentioned by Geballe & Van der 
Veen (1990) suggesting that it is an evolved F-type super- 
giant displaying carbon-rich circumstellar dust. Its 21/im 
feature was discovered by Kwok et al. (1995). An ex- 
tensive abundance analysis was previously carried out by 
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Table 5. Chemical analysis of the programme stars. For every ion we list the solar abundance and for every star the number of lines used (N), t 
equivalent width (W\), the abundance ratio relative to iron ([el/Fe]) and the internal scatter (a ), if more than one line is used. For the solar iron at 
we used the meteoric iron abundance of 7.51. For the solar C, N and O abundances we adopted resp. 8.57, 7.99 and 8.86 (C: Biemont et al. 
Hibbert et al. 1991, O: Biemont et al. 1991); the other solar abundances are taken from Grevesse (1989). Note that abundances based on one lin 
be treated with caution. 
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Fig. 4. The abundances of our six programme stars relative to iron [el/Fe] 



For the model parameters, we found Tcff=6500K, 
log5=1.0 and ^t=3.5 kms~^. The model parameters agree 
with those found by Reddy (6500, 0.5, 5.0), but there are 
some differences in the abundances. In general, for species 
displaying a large difference, our analysis was based upon 
significantly more lines, which should imply that our re- 
sults are intrinsically more reliable. 

Besides the huge carbon enrichment also oxygen is 
moderately enhanced ([O/Fe]=-|-0.6) but this value is 
based only on two blended lines. 

Reddy et al. also found a significant overabundance 
of lithium and we confirm the presence of the Li line at 
A670.7nm in our spectra as well. 

The number of lines of the s-process elements is im- 
pressive. A close look at the spectrum reveals that it is 
even completely dominated by lines of s-process elements. 
[s/Fe] (Table |) is 2.2, making IRAS05341 the most s- 
process enhanced post-AGB star known so far. Note the 
even-odd alternation of the absolute abundances (Fig. ^ 
and the good internal consistency {a < 0.20). There is 
significant difference in the [el/Fe] ratios of our analysis 
and the one of Reddy et al. (1997) reaching up to -1-0.8 
dex. Since we take much more and much weaker lines into 
account reaching a good internal accuracy, and use the 



5.3. IRAS 07134+1005 = HD 56126 

With its 8.3 visual magnitude, this is the brightest star 
in our sample. Considerable efforts have been performed 
already to reveal the properties of the circumstellar shell 
of this well-known post-AGB star. i^CQ, i3cO, HON and 
HCO+ (Bujarrabal et al. 1992, Omont et al. 1993) have 
been detected in millimeter and radio line emission. The 
21//m feature was discovered by Kwok et al. (1989) in its 
IRAS Low Resolution Spectrum. Optical spectra revealed 
the presence of circumstellar C2, ^^CN and -^^CN bands 
(Bakker et al. 1996; Bakker & Lambert 1998a) and i^ci^C 
and i^ci^O bands (Bakker & Lambert 1998b) , which led 
to an estimate for the i^c/iSC ratio : i^c/isc^ 72 ± 
26. Mid-infrared (8-21/im) imaging shows no spherical but 
possibly an axial symmetry of the dust emission (Meixner 
et al. 1997; Dayal et al. 1998). 

Concerning the photosphere, there is general agreement 
that IRAS07134 is a pulsating star, but different pulsa- 
tion periods have been published (e.g. Lebre et al. 1996). 
Limited abundance analyses have already been performed 
by Parthasarathy et al. (1992) and by Klochkova (1995), 
showing chemical evidence for the post-AGB status of this 
object. However, the high quality NTT-spectra for this 
star permitted us to elaborate a much more detailed abun- 
dance analysis. We used 342 lines of 27 different ions. 
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However, the latter author found a much larger abundance 
if deduced from lines of neutral atoms than from lines of 
singly ionized atoms of the same element. This was at- 
tributed to unknown non-LTE effects. We do not confirm 
this result since careful analysis of our much higher S/N 
spectra does not reveal the presence of most neutral lines, 
while we identified some lines differently. This illustrates 
that it is often very difhcult to assign uncertain LTE abun- 
dance results of an individual object solely to non-LTE ef- 
fects. Also the effect of other uncertainties should be care- 
fully studied and a more systematic approach using data 
on several objects should be initiated to evaluate possible 
non-LTE deviations of specific lines. 

The ionization equilibrium of Ni is also fulfilled but for 
Cr we find a A(log A(Crll) - logA(Crl)) = 0.34. A lower 
model gravity can solve this inconsistency for Cr but it 
may also indicate non-LTE occupation numbers for the 
different ions of Cr. We chose not to change our model 
parameters in order not to compromise significant com- 
parison with the other stars in our sample. 

A problem occurred when calculating the Fe abun- 
dance from lines of neutral iron with a very low excita- 
tion potential: the lines with a very small lower excitation 
potential gave systematically lower abundances. This is 
a strong indication that non-LTE over-excitation occurs. 
We therefore discarded all FeLlines with an excitation po- 
tential X < 1.25eV in our abundance determinations, also 
for all other programme stars. 

The large line-to-line scatter in the abundance of Y and 
Zr is solely caused by only three fines of Y (A578.169nm, 
A550.990nm and A528.982nm), and two lines of Zr 
(A448.545nm and A449.546nm). These lines do not show 
any systematic trend with another atomic parameter 
(wavelength, excitation potential, equivalent width), so 
the anomalies are probably caused by inaccurate log{gf)- 
values. The only line in common with the other stars is 
the Y II- line at 528.982 nm in IRAS22223. Also for this 
star, the abundance inferred from this line is significantly 
lower than the abundances inferred from the other lines. 

We also detected several lines of singly ionised 
Hafnium (Z=72) in the blue part of the spectrum of 
IRAS07134, but only two of them turned out to be un- 
blended (A409.316nm, Wx = 83 mA and A466.414nm, 
Wx = 19 mA). Unfortunately, the Hf-abundance derived 
from these two lines differs by 0.4 dex (logA(Hf) = 1.72 
and 2.09 respectively). These abundances are both very 
high if compared with the calculations of Malaney (1987). 
We expect logA(Hf) « 1.1 for a neutron exposure of 
T = 0.3 mb^^ (neutron exposure: see Sect. 0). No other 
possible identification could, however, be found. An in- 
terpretation of this result is difficult, not only because of 
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Fig. 5. The HflLhne (A409.316nm) in the blue spectrum 
of IRAS07134 and IRAS19500. The derived hafnium abun- 
dance is quite high, but blends or a wrong identification 
is almost certainly excluded. 

5.4. IRAS 19500- 1709 ^ HD 187885 

This post-AGB candidate shows only a weak 21/im feature 
(Justtanont et al. 1996) which is probably an effect of the 
higher temperature of this object. We already pubfished 
an abundance analysis of this object (Van Winckel et al. 
1996a) based on ESO La Silla spectra. The latter spectra 
show Hcl lines at 447.08 nm and 587.6 nm pointing to a 
hot source. Indeed, demanding excitation equilibrium of 
the Fel-lines yielded a temperature Toff=8000K. 

In this study we present a new analysis (Table ph based 
on UES-spectra with a higher S/N and a broader spectral 
coverage. The analysis of the iron lines yielded the same 
atmospheric parameters, only the log(iyA/A)-abundance 
diagram points to a larger microturbulent velocity ^t — 
6km s^^. The lines of the the s-process elements Y and Zr 
favour a slightly lower value of £,t < 5km s""'^. We again 
keep the atmospheric parameters obtained using the iron 
lines. A change in the microturbulent velocity of ^Ikms^^ 
will only infiuence the abundance deduced from the largest 
lines {Wx > 120mA) with a typical value of 0.1 dex. 

Comparing this analysis with our previous one, there 
are some quantitative differences, but our conclusions for 
this star remain unchanged. The s-process overabundance 
is now much better established, with more lines for the 
Is-elements and new abundances for the hs-elements La, 
Ce, Nd and Eu. 

We report also for this star the detection of a line 
of singly ionised Hafnium at A409.316nm. An equivalent 
width Wx — 12 mA gives an abundance of logA(Hf) — 
1.53 which is again quite high if compared with the cal- 
culations of Malaney (1987). We expect log^(Hf) w 1.0 
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Table 6. S"^ values for different neutron exposures tq 
(N„ — lO^cm"'^). The stars are listed with increasing 
[hs/ls]-value. 
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should have an equivalent width of Wx = 9 mA. Inac- 
curate log(g/)-values can be the cause of this 
inconsistency. Other possibly detectable Hf II lines are all 
heavily blended. 

5.5. IRAS 22223+^327 

The chemical composition of this star was recently dis- 
cussed by Decin et al. (1998). We list the individual abun- 
dances here again since minor differences are obtained due 
to an upgrate of the line-list and a homogenization of the 
atomic data. 



5.6. IRAS 23304+6147 

Kwok et al. (1989) classified this object as a proto-pla- 
netary nebula and discovered its 21/im feature. The cir- 
cumstellar dust turned out to be C-rich as confirmed by 
several authors: detection of C2, C3 (Hrivnak 1995), CO, 
HCN (Omont et al. 1993), CN (Bakker et al. 1997) and 
non-detection of OH (Likkel 1989). Hrivnak et al. (1999) 
found the object to be resolved in V with an extent of 
about 2 arcsec. Reddy & Parthasarathy (1996) derived 
model parameters for this object by fitting a Kurucz model 
atmosphere through the observed spectral energy distri- 
bution. This technique yielded Teff=5000K and log 3=1.0. 
Comparing our spectra with the other stars in our sam- 
ple, we concluded that this temperature was probably too 
low : the high excitation 0-triplet around A615.6nm is 
clearly detected (see Fig. |^). Indeed, also our spectroscopic 
determination of the temperature using the iron-lines re- 
sulted in a Teff=6750K. The abundance analysis makes 
IRAS23304 join the group of post-AGB stars that clearly 
display chemical evidence of their post 3rd dredge- up char- 
acter. 

To our knowledge, this is the first detailed photo- 
spheric abundance study of this star. The internal scat- 
ter of the four 0-lines is small indicating that the very 
high C/0 ratio of 2.9 is probably real. For the metallic- 
ity, we find the value [Fe/H]=— 0.8. The other iron peak 
elements do follow this deficiency. Again, the s-process 
abundances are quite impressive for this star. Note that 
the internal scatter of the different abundance determina- 
tions are somewhat higher for this star than for the other 
programme stars as a consequence of the lower S/N of our 
spectra (m„ = 13.1). 

6. Neutron exposure of 21/xm stars 

6.1. Goodness-of-fit 

A theoretical parameter describing the neutron irradia- 
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modelled by an exponential irradiation parameterized by 
a mean neutron exposure tq (e.g. Ulrich 1973). In most 
calculations, the neutron density and the fraction of the 
inter-shell material that remains exposed to neutrons are 
adopted as constants. The efhciency of s-processing is then 
indicated by one parameter tq. 

The neutron exposure tq is intimately related to the 
[hs/ls] index: a large tq yields a large [hs/ls]. In order to 
estimate tq quantitatively, we used the theoretical models 
of Malaney (1987) and the 'goodness of fit' procedure as 
defined by Cowley & Downs (1980). A description of this 
method can also be found in Vanture (1992) and Smith et 
al. (1996). In this method, the quantity S"^ is defined as 



1 



S' = i- 



N 



io,-iM, + s)y 



N ^^ 

i=l 



where N is the number of s-process elements involved in 
the comparison, Oi the observed abundance of element 
i, Mi the predicted model abundance as tabulated by 
Malaney, S the average offset between the observed and 
the predicted s-process abundances {S — ^/NJ2i {Oi — 
Mi)) and (Tj is the uncertainty on the obtained individ- 
ual abundance Oi. Oi and ai are calculated as described 
in Decin et al. (1998) so on top of the internal accuracy 
a fixed value of 0.3 was quadratically added to account 
for the uncertainties on the model atmosphere parame- 
ters and eventual systematic offsets in the log(gf) values. 
The values for S^ for each model and each star are listed 
in Table o. A graphical presentation of the best fit for 
two stars (IRAS07134 and IRAS05341) can be found in 
Fig. 0. Note that the S^ are small which indicates a good 
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Fig. 7. The correlation between the total enrichment of 
s-process elements ([s/Fe] which is the mean of [Y/Fe], 
[Zr/Fe], [Ba/Fe], [La/Fe], [Nd/Fe] and [Sm/Fe]) and the 
[hs/ls] index. The straight line gives the simple least- 
squares fit. 



Fig. 6. Comparison of the observed abundances with 
the calculations of Malaney (1987) for IRAS07134 and 
IRAS05341. The observed abundances were corrected for 
the initial abundances (A^) assuming [s/Fe].i=0 (except 
for Eu). As evidenced by the 'goodness of fit' procedure, 
the s-process elemental distribution of IRAS07134 points 
to a rather moderate neutron exposure, whereas the dis- 
tribution of IRAS05341 is best fitted by a large neutron 
exposure tq = 0.8 mb^ . 



S'^ and so the neutron exposure tq is well defined. For stars 
with a higher [hs/ls], a larger range of values for tq are 
possible: a high neutron exposure will only raise the overall 
level of enhancement without significantly changing the 
distribution of the individual abundances (see also Fig. 
P). Therefore, the s-process distribution of IRAS05341 fits 
well for every tq > 0.5 mb^ . 



6.2. Total s-process enrichment of 21^m stars. 

In Fig. ^ the [hs/ls] index from the 21/xm stars is shown 
in comparison with the overall s-process enrichment as 
parameterised by the [s/Fe]-index. More objects should 
clearly be studied, but the total neutron exposure is sur- 
prisingly strongly correlated with the total enrichment of 
s-process elements in the 21/im stars, with more enriched 
objects showing also a stronger total neutron irradiation! 
A simple linear least-squares fit gives a high correlation 
coefficient of -t-0.96 with the fit 



1995), the [hs/ls] index will not change with increasing 
dredge-up of enriched material in the stellar envelope, 
making the [hs/ls] index independent of the dilution fac- 
tor of the enriched material by the initial stellar envelope 
material. The [hs/ls] correlation with [s/Fe] therefore re- 
flects a different and more efficient internal nucleosynthe- 
sis, with increasing dredge-up efficiency ! In AGB evolu- 
tionary models, the inclusion of protons into the intershell 
will create a primary -'^^C pocket and induce the neutron 
production. Fig. indicates therefore that, with increas- 
ing dredge-up efficiency, also the dredge-in of protons will 
increase, giving rise to a higher neutron irradiation. 

6.3. Metallicity. 

Although our sample at this point only consists of 6 ob- 
jects, a considerable spread in metallicity ranging from 
—0.3 to —1.0 is observed. Since the ^^C(a,n)^^0 neutron 
source in AGB stars is most likely of primary origin the 
[hs/ls] index should increase with decreasing metallicity 
due to the larger amount of neutrons per iron seed nucleus 
(e.g. Clayton 1988). This assumes, however, that the dif- 
fusion of protons into the C-rich layer itself is independent 
of metallicity or on parameters linked to the metallicity 
(like temperature profile, luminosity etc.). A higher neu- 
tron exposure is therefore expected for low metallicity ob- 
jects assuming that the thermal pulse phenomenon itself 
is independent of metallicity. 

In Fig. pi we display the [hs/ls]-index in function of the 
metallicity. The 21/im stars clearly cover a wide range in 
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bert (1985,1986,1990) and Smith et al. (1987) : since the 
spectra of oxygen rich AGB stars are dominated by strong 
TiO-bands, only a few spectral windows can be used to in- 
fer abundances from atomic lines. The regions in between 
740 and 758 nm, and 998 to 1010 nm, were used together 
with additional IR spectra around 1.6, 2.2 and 4.0 /im 
to deduce abundances of CNO isotopes. The situation is 
even worse for C stars where the C2 and CN bands domi- 
nate the optical and IR spectrum. To our knowledge only 
Utsumi (1985 and references therein) made a quantitative 
analysis of the s-process enhancement in C stars using 
two small spectral windows with minimal molecular ab- 
sorption (475-490 nm, and 440-450 nm). SC stars, on the 
other hand, are AGB stars with a C/0 ratio close to 1.0 
making the oxide bands and bands from carbon molecules 
very weak or absent so that many atomic lines become de- 
tectable in the spectrum. We use the results for seven SC 
stars from a recent paper by Abia & Wallerstein (1998). 

In Fig. we compare the observed overabundances of 
the light s-process elements ([Is/Fe]) with the ones ob- 
served in intrinsic AGB stars. Since there is strong ob- 
servational evidence that the 21/im stars are post-carbon 
stars, it is no surprise that the overabundances of the s- 
process elements in 21/im stars are indeed large in com- 
parison with MS, intrinsic S and SC stars and are com- 
parable with and on average even slightly higher than ob- 
served in C stars. We therefore can confirm earlier results 
that the increasing C/0 ratio in the MS-S-SC-C sequence 
is also reflected in an increasing enhancement of s-process 
elements (e.g. Smith & Lambert 1990; Abia & Wallerstein 
1998). 

1.2. Extrinsic enriched objects 

To investigate the metallicity dependence of the integrated 
neutron exposure for intrinsic AGB stars is not easy since 
the metallicity of field carbon stars is very difficult to de- 
termine and the spread in metallicities observed in Galac- 
tic M-MS-S-SC stars is not large (e.g. Smith & Lambert 
1990). Plez et al. (1993) have shown, however, that seven 
intrinsic AGB stars in the metal poor ([Fe/H] = —0.5) 
Small Magellanic Cloud do indeed show a very high neu- 
tron exposure and they suggest that it is not intrinsic to 
the SMC but due to the low initial metallicity. 

The metallicity spread observed in the extrinsic en- 
riched objects is much larger and there is growing ob- 
servational evidence that the [hs/ls] index is indeed anti- 
correlated with metallicity (e.g. Smith 1999): Busso et 
al. (1995) confirmed earlier suggestions by Luck & Bond 
(1991) in a compilation of literature values of the [hs/ls] in- 
dex in Ba stars and MS-S stars and quantified tentatively 
the metallicity dependence of the neutron irradiation as 
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Fig. 8. The [hs/ls] index of the analysed 21/im stars as 
a function of the metallicity determined by the Fe abun- 
dance ([Fe/H]. The line gives the simple least-squares fit 
[hs/ls] = -0.93 [Fe/H] - 0.75. 



is -0.67). The most metal deficient object IRAS07134 has 
seem to have undergone a much smaller neutron irradia- 
tion than the less metal deficient IRAS05341. Clearly, also 
other stellar parameters determine the internal AGB nu- 
cleosynthesis during the AGB evolution. 

Note that the two stars significantly below the least- 
square line are the ones for which we detected the very 
heavy element Hf (Z=72). If indeed these high abundances 
turn out to be true, this would indicate that the s-process 
is able to produce elements heavier than the Ba peak and 
thus probably up to the lead-peak! If true, the low [hs/ls]- 
index of the two stars might be due to the production of 
very heavy elements. It is clear that high S/N blue spectra 
of all the objects are needed for a more detailed discussion. 

Moreover, several (probably single) post-AGB objects 
are known which do not show s-process enrichment nor 
enhanced C/0 ratios. These objects include HD 161796 
([Fe/H]=-0.3; Luck et al, 1990), HD 133656 ([Fe/H]== 
-1.0; Van Winckel et al. 1996b) and SAO 239853 ([Fe/H] = 
—0.8, Van Winckel 1997) which show very similar metal- 
licities and IR-excesses. Unfortunately, like the 21/.tm ob- 
jects they are too far for reliable HIPPARCOS parallaxes 
so the luminosity differences between the enriched and 
non-enriched objects is an as yet unknown but probable 
key parameter in this discussion. 



7. Comparison with intrinsic and extrinsic 
enriched stars. 

7.1. Intrinsic enriched objects: AGB stars 
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Fig. 9. Comparison between the overabundances of light 
s-process elements observed in the 21/xm stars with those 
from intrinsic AGB stars. 'Is' is the mean of the light s- 
process elements Y and Zr, while 'hs' is the mean of the 
heavy s-process elements Ba, La, Nd and Sm. The refer- 
ences for the observations of the intrinsic AGB stars can 
be found in Abia & Wallerstein (1998). The full lines are 
the theoretical predictions for different neutron exposures 
paranieterised by tq (in mbarn^^) arc from Busso et al. 
(1995). 

larging the metallicity spread to lower metallicities. North 
et al. (1994) also found that the [hs/ls] index was corre- 
lated with the metallicity in a sample of CH-subgiants. 
Finally also the low metallicity yellow symbiotics, who re- 
cently were included in the Ba-star family, show s-process 
distributions characterised by a relatively high neutron ir- 
radiation (Smith et al. 1996, 1997; Pereiraet al. 1998). The 
results obtained by Luck & Bond (1991) included, how- 
ever, 4 metal deficient objects with small overabundances 
showing a very weak neutron irradiation. They dubbed 
these stars "the metal-deficient Ba stars" . 

Some caution should be expressed by generalising these 
results to single star AGB evolution : especially for the 
short orbital period binaries there is strong evidence that 
the binarity affects the mass-loss and thus the AGB evolu- 
tion itself of these binaries. The amount of thermal pulses 
and the characteristics of the s-process synthesis might 
therefore also be affected by the orbital elements (e.g. 
Jorissen & Mayor 1988). 

In Fig. no we compare the [hs/ls] index in function of 
the metallicity with the values found in the literature for 
other s-process enriched objects. This compilation is diffi- 
cult since the definition of the hs and Is indices is author- 
dependent. We scaled the abundances to the definition 
given in section 4.4. For undetermined abundances of the 
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Fig. 10. The [hs/ls] index as a function of the metallic- 
ity for a sample of s-process enriched objects. Asterisks 
are the 21/im stars; triangles the metal deficient Barium 
stars from Luck & Bond 1991; squares the CH-giants from 
Vanture (1992); circles the Ba-stars from the compilation 
made by Busso et al. (1995); plus-signs the CH-subgiants 
from North et al. (1994); crosses the yellow symbiotic stars 
from Smith et al. (1996,1997) and Pereira et al. (1998); di- 
amonds the SC stars form Abia & Wallerstein (1998). 

in total 7 elements (Sr, Y, Zr, Ba, La, Nd and Sm) were 
available. 

Fig. |l^ shows that there is a large spread in the in- 
tegrated neutron irradiation, but that there is indeed a 
global increase with decreasing metallicities from [Fe/H] 
~ —0.5 downwards as already noted by several authors 
(e.g. Smith, 1999). Three out of four "metal-deficient Ba 
stars" show a very low [hs/ls]-index (Luck & Bond, 1991). 
These stars certainly are worth a re-analysis based on 
higher resolution, higher signal-to-noise spectra since the 
internal errors of the analysis based on photographic in- 
termediate resolution (0.2 A) spectra are rather large and 
the observed overabundances small (Luck & Bond, 1991). 
Some CH-stars show extremely high [hs/ls] -indices for 
their metallicity but also some Ba stars like the very Nd- 
rich HD 35155 (not in the figure) show a-typical values. 
The 21/im post-AGB objects follow the same trend as 
the CH-subgiants of the same metallicity and show on av- 
erage a somewhat lower [hs/ls]-index compared to the Ba 
giants. 

We can conclude that the metallicity dependence of 
the neutron exposure has a large intrinsic spread which 
confirms that certainly also other fundamental parame- 
ters determine the internal nucleosynthesis and dredge-up 
phenomena during the AGB evolution. 
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With their low metaUicity ([Fc/H] ranging from —0.3 to 
— 1.0), high C/0 ratio (ranging from 1.0 to 2.9) and espe- 
cially their high photospheric abundance of s-process el- 
ements, we confirm that they are post-3rd-dredge-up ob- 
jects. Since also the circumstellar environment is C-rich 
and no orbital motion is observed we can classify the 21^m 
stars as single post-carbon stars of low initial mass. 

The G-F spectral type makes that a wide variety of 
s-process elements can be detected by their atomic lines. 
This, together with their significant spread in metaUicity, 
make that the 21/xm stars are ideal objects to constrain 
the AGB nucleosynthesis and dredge-up models. The most 
interesting constraint in this context comes from the pos- 
sibility to accurately determine a wide variety of elemental 
abundances of these intrisically enriched objects. Contrary 
to AGB stars themselves, no photospheric molecular ab- 
sorption is observed, so isotopic ratios of key elements are 
in the 21/im stars unfortunately difficult to determine. 
The molecular bands originating in the circumstellar en- 
vironment can, however, be used for some isotopes like 
i^C and "C (e.g. Bakker & Lambert 1998a). 

We can summarise the main conclusions with respect 
to the s-process nucleosynthesis displayed by 21/im stars 
as follows : 

— The 21/im stars display a wealth of atomic lines 
from s-process elements. The photospheric spectrum 
of stars like IRAS05341 is even completely dominated 
by atomic lines of heavy elements. The observed over- 
abundances of s-process elements in the analysed 21/im 
stars are large and range from [Is/Fe] = +1.3 to +1.9 
which make the 21/im stars among the most s-process 
enriched objects known. 

~ The 21/xm stars display a wide variety of integrated 
neutron irradiation as parameterised by the observa- 
tional [hs/ls]-index, corresponding to a tq neutron ex- 
posure from To = 0.2 to 0.8 mbarn"^. 

— The [hs/ls] index displayed by the 21/im stars is 
strongly correlated (classical correlation coefficient of 
+0.96) with the total enrichment of the s-process el- 
ements as parameterised by the [s/Fe]-index, in the 
sense that more enriched objects also display a higher 
integrated neutron irradiation. Since in carbon stars 
the asymptotic values of the s-process distribution is 
probably reached, this means that the dredge-up effi- 
ciency is strongly linked with the neutron production 
in the intershell ! 

— The anti correlation of the [hs/ls] index with the metal- 
licity of the 21/im objects is less well determined 
(classical correlation coefficient —0.67) and contains 
a large scatter. This scatter is certainly intrinsic and 
confirms that also other fundamental parameters de- 
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— Compared with other s-process enriched objects of ex- 
trinsic nature, the neutron exposure displayed by the 
21/im stars fall roughly in the same range as the CH 
subgiants and is on average lower than it is in Ba- 
giants of the same metaUicity. Note that the metaUic- 
ity dependence of the neutron exposure displayed by 
intrinsic and extrinsic objects contains also a large in- 
trinsic scatter. Especially the iron deficient Ba stars 
as defined by Luck & Bond (1991) should be analysed 
in great detail since they display a very low neutron 
irradiation for their low metaUicity. 

We can conclude that a homogeneous detailed study of 
the 21/im stars forms a very useful complement to other 
types of s-process enriched objects (intrinsic AGB stars 
and extrinsic enriched objects) used in the observational 
study of the theoretically less well understood s-process 
nucleosynthesis, mass-loss and dredge-up phenomena oc- 
curring on the AGB. 
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